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Abstract 

We use digital holographic microscopy to measure the three-dimensional (3D) diffusion tensor 
D for clusters of colloidal spheres, including both dimers and triangular trimers of spheres, for 
which no analytical calculation of the diffusion tensor exists. We measure all nonzero elements 
of the diffusion tensor to a precision of 1% or better, and we detect anisotropic rotational and 
translational diffusion arising from the asymmetries of the clusters. In the case of the three- 
particle triangular cluster, we also detect a small but statistically significant difference in the 
rotational diffusion about the two in-plane axes. We attribute this difference to weak breaking of 
threefold rotational symmetry due to a small amount of particle polydispersity. Our experimental 
measurements agree well with theoretical predictions and show how diffusion constants can be 
measured under conditions relevant to colloidal self-assembly, where theoretical and even numerical 
prediction is difficult. 
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Diffusion plays a critical role in the dynamics, self-assembly, and rheology of complex 
fluids. In systems such as colloidal suspensions, which typically have short-ranged interaction 
potentials, diffusion can in fact play a larger role than energy barriers in setting transition 
rates However, the diffusion of geometrically anisotropic particles, a common class of 
colloidal suspension that can also arise as intermediates in the self-assembly of spherical 
particles, can be difficult to predict. Theoretically determining friction factors for these 
particles requires analytically solving Stokes' equation, which is only possible for highly 
symmetric particles such as ellipsoids [2] or sphere dimers [3] in unbounded fluids. Numerical 
methods such as bead modeling [4j or finite-element methods [SHZ] require approximating the 
shape of the particles or the hydrodynamic interactions. These methods are difficult to apply 
when asymmetric particles diffuse near rigid boundaries or other particles, two situations 
that are relevant to colloidal self-assembly and dynamics in general. Thus experimental 
measurements of diffusion tensors are crucial. 

In particular, measurements on single particles rather than ensembles are necessary. 
Anisotropic particles show multiple diffusion timescales that are difficult to resolve through 
bulk techniques such as depolarized dynamic light scattering |8]. Existing techniques for 
measuring single-particle diffusion are also difficult to apply in 3D. Video microscopy, for ex- 
ample, has been used to measure two-dimensional (2D) diffusion of colloidal ellipsoids PfTU] 
and planar sphere clusters [TT] but yields limited information about out-of-plane motions 
[T2HTi] . Confocal microscopy can be used to study the 3D Brownian dynamics of a range 
of geometrically anisotropic particles [151 HB] , but can take on the order of 1 s to acquire a 
3D stack. This can make it challenging to probe timescales comparable to particle diffusion 
times or to study rare processes such as the early stages of self-assembly [T7] . 

In this Letter, we present precise measurements of 3D diffusion tensors of individual 
colloidal clusters using a fast 3D imaging technique, in-line digital holographic microscopy 
(DHM). In DHM one records a 2D hologram generated by interference between light scat- 
tered from colloidal particles and the undiffracted, transmitted beam (Fig. [ij^a)). The 2D 
holograms can be acquired much more rapidly than 3D confocal stacks. By fitting models 
based on electromagnetic scattering solutions [T8H20], we recover the 3D dynamics of dimers 
and triangular trimers of colloidal spheres. We resolve all the translational and rotational 
components of the diffusion tensor and detect the asymmetries of these clusters. We also 
show that a small amount of polydispersity leads to symmetry breaking in the rotational 
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FIG. 1. (Color online) (a) Schematic overview of in-line holographic microscope used in the 
experiments, (b) Axis trajectories for a diffusing triangular cluster of 1 /xm diameter polystyrene 
spheres (inset: orientation of cluster axes). Renderings show the trajectories of axes 1 and 3 on 
a unit sphere. Blue (dark) markers indicate the start of the 40 s trajectory and Orange (light) 
markers the end. (c) Laboratory-frame center-of-mass trajectories of the same triangular cluster 
over the same time interval as in (b). Inset: amplitude reconstructions [21] of the holograms, 
showing the structure of the measured clusters. Scale bar, 2 /xm. 

diffusion components. Our measurements agree well with numerical predictions and, more 
generally, show how diffusion tensors can be measured in experimental systems relevant to 
self-assembly, where theoretical predictions are challenging. 

The diffusion tensor D quantifies the translational and rotational Brownian diffusion of 
an arbitrary rigid colloidal particle. In general, six generalized coordinates (three positions 
and three orientation angles) are needed to describe the position and orientation of a rigid 
body. Although D is rigorously defined by generalizing Pick's law to an abstract ensemble 
of particles diffusing in this 6-dimensional configuration space, D also describes correlations 
between displacements of the for short lag times r [221 123] [21] : 
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D may be partitioned into the following 3x3 blocks: 

D = ( • (2) 

D** describes translational diffusion, D^^ rotational diffusion, and D*^ translation-rotation 
coupling. The generalized Stokes-Einstein equation relates D to the friction tensor ^ de- 
scribing the hydrodynamic Stokes drag forces and torques on a moving particle: D = 
ksT^'^ [221 [23]. Predicting D thus requires a solution for the Stokes flow around a particle. 
Dimers are one of the few non-spherical shapes for which analytical solutions exist; none 
exist for trimers. 

We make dimer and triangular trimer clusters through limited aggregation of sulfate 
polystyrene spheres (Invitrogen) |25], 1.3-/im diameter for dimers and l-/im diameter for 
trimers. We transfer these particles into a 250 mM NaCl solution to screen the charge of the 
stabilizing sulfate groups and start the aggregation, then we decrease the ionic strength by 
quenching with deionized water (Millipore) after 1 minute to arrest the aggregation. We then 
suspend the resulting mixture of single particles, dimers, and larger clusters in a density- 
matched solvent consisting of 50% v/v D2O and 50% v/v H2O with a salt concentration 
of 1 mM. We load the particles into sample cells made from glass slides, coverslips, and 
76- /im thick Mylar spacers [191 ED]. After finding a cluster with the desired morphology 
using bright field microscopy, we ensure that it is at least 30 /im away from sample cell walls 
or other particles so that the cluster diffusion is unhindered by boundaries [26]. We record 
holograms using an instrument previously described in the literature [TTJ |20] at a frame rate 
of 25 frames per second and a temperature T = 296 ± 2 K. 

We then analyze the measured holograms to obtain 3D trajectories of the clusters. We 
obtain particle sizes, refractive indices, center-of-mass 3D positions, and three orientation 
angles (two for dimers) by fitting an exact scattering solution to Maxwell's equations to 
each recorded hologram [m [27|[28]. Figs. |l|(b) and (c) show some of the 3D data we 
obtain for a trimer. We measure the components of the translational block D** by directly 
applying Eq. [T| where the relevant Ag' are relative to a coordinate system rigidly fixed to the 
cluster [ini [22l [23] . The correlation functions needed to measure the diagonal components 
of D** are cluster- frame mean-squared displacements (MSDs). To measure the components 
of D'^'', we examine the dynamics of the axis vectors Uj fixed to a cluster. The tips of 
the Uj diffuse along the surface of a unit sphere, as illustrated in Fig. [l|^b). We compute 
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autocorrelations of the Uj, which are related to -0^,1, the 3 diagonal components of D'"'', as 
follows 



{ui(t) ■ Uj(t + r)) = exp 




(3) 



For the axisymmetric dimer, we consider a related quantity, the MSD of the axis unit vector 
u [ISlIISlEniEI]: 

(Au2(r)) = 2(l-(u(t).u(t + r))). (4) 

We find good agreement between D for a dimer, measured to 0.5% precision from a 
time series of 22,000 holograms, and analytical and numerical predictions. The measured 
axis MSD, along with a best fit to Eq. |4| and the cluster-frame MSDs are shown in Fig. |2} 
Dimers of two identical spheres have three mirror planes and an axis of continuous rotational 
symmetry, so in the coordinate system shown in the inset of Fig. [2](b), D is diagonal. But 
due to the breaking of spherical symmetry, D has four rather than two unique elements: 
the translational diffusion constants and D±, and the rotational diffusion constants Dr^\\ 
and -Dr,±; where the subscripts denote motions parallel to and perpendicular to the dimer 
axis [Tnil22]- Because we cannot observe rotations about the dimer axis, we can only measure 
Dr^±. Moreover, we can only measure a combined MSD, {Ax±{t)'^), representing translation 
along the two degenerate perpendicular axes. We then extract and D± from linear fits 
to the MSDs. (Axy(r)) has a slope of 2D\\ in accordance with Eq. jT| and {Ax']_{t)) has a 
slope of 4:D±. Our measurement of D\\/D±, which is a universal constant for any Brownian 
dimer, agrees well with predictions from the shell modeling code hydrosub [32j and the 
exact Stokes solution of Nir & Acrivos [3J (Table [l]). 

The individual elements of D also agree with calculations, once we account for the sphere 
radius a and the solvent viscosity rj. We determine an effective sphere radius Oefr from the 
measured ratios D^\/Dr^± and D±/Dr^±, which depend only on a [3]. Then we determine the 
best-fit t]es from the measured elements of D. We find an effective acs = 709 nm, which is 
larger than the optical radius Oopt = 650 nm obtained from DHM. The larger effective ra- 
dius is consistent with typical dynamic light scattering measurements of the size of colloidal 
spheres, which show enhanced hydrodynamic radii due to charge or hairy surface layers on 
the particles [MIES]. The best-fit viscosity is r^ofi = 1.159 mPas, consistent with measure- 
ments of the diffusion constant of a freely diffusing sphere in the same sample, which acts 
as an in situ thermometer, combined with bulk viscosity measurements using a Cannon- 




r(s) 

FIG. 2. (a) Axis MSD for dimer of 1.3-^m diameter spheres. Open symbols are measurements; 
solid line is the best-fit to Eq. |4| where (u(t) • u{t + r)) = exp(— 2L>j. _|_r). (b) Cluster-frame MSD 
computed for the same dimer as in (a). Open symbols are measurements; solid lines are linear fits. 
Triangle shows slope of 1. Error bars are calculated using a block decorrelation technique [33]; 
they are comparable in size to the plotting symbols or smaller. Inset: orientation of parallel (||) 
and perpendicular (_L) axes. 

Manning capillary viscometer; these together give a solvent viscosity of 1.19 ± 0.04 mPas. 
We use this procedure because the solvent viscosity has a strong temperature dependence 
[28] . The elements of D computed with these effective parameters agree with our measure- 
ments to better than 1% (Table [l]). We also note that the hydrosub prediction for -Dr,± 
differs from the analytical prediction by about 1%, which is consistent with prior studies [1]. 
The agreement between our measurements and the analytical prediction suggests that our 
measurement accuracy is at least comparable to, if not better than, that of hydrosub. 

Measurements on trimers reveal anisotropic translational and rotational diffusion. Trimers 
of identical particles have two mirror planes, making D** and D'^'^ diagonal [221 EB] in the 
coordinate system shown in the inset of Fig. [1]d. We denote the six diagonal elements 
as -Dt,i, -Dt,2, -Df,3, -Dr,i, -Dr,2, and -0^,3. In contrast to dimers, trimers lack axisymmetry, 
allowing us to observe rotations about all three axes and measure all six elements. In Fig. |3j 
we show the axis autocorrelations (uj ■ Ui(t + r)) computed from 20,000 holograms, as well 
as best fits to exponential decays. The autocorrelation of axis 3 decays more rapidly than 
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TABLE I. Measured diffusion tensor elements for dimer in Fig.[2| along with analytical calculations 

from an exact Stokes solution [3j and numerical calculations from hydrosub [32]. Calculations 
use a best-fit particle radius Ocff = 709 nm and solvent viscosity rj^s = 1.159 mPas. 

Experiment Exact hydrosub 

Dr,^_ (s-i) 0.1034 ± 0.0006 0.1034 0.104 

D|| (xl0"^3 m^s-^) 2.015 ± 0.012 2.010 2.02 

-D± (x 10-13 i^2g-i) 1.785 ± 0.007 1.790 1.80 

Du/Dj_ 1.129 ± 0.011 1.123 1.12 




FIG. 3. Cluster axis autocorrelations (uj(t) • Uj(t + r)) for a trimer of l-/im diameter spheres, 
showing anisotropic rotational diffusion. Open symbols are experimental measurements; error bars 
are comparable to or smaller than symbols. Solid lines are fits to exponential decays. Inset shows 
cluster axis orientation. 

the autocorrelations of axes 1 and 2, in agreement with expectations: as shown in Eq. |3j 
(u3(t) ■ u^it + t)) depends on D,. i and -Dr,2, both of which should be larger than 1)^,3 due 
to hydrodynamics. The elements of the diffusion tensor that we extract from this data are 
shown in Table [Tij The difference between 3 and both Dr^i and Dr^2 is much larger than 
the experimental uncertainty, showing clear evidence for anisotropic rotational diffusion. 
The translational diffusion we observe is similarly anisotropic (Fig. |4]and Table [ll|) . 

Interestingly, although our measurements of the dimensionless ratios Dr^i/Dr^s and 
Dt,i/Dt;i agree well with the HYDROSUB predictions, we observe small but statistically 
significant differences between the elements of D corresponding to the two in-plane axes 1 
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FIG. 4. Body-frame MSDs for the same trimer in Fig. |3| See inset in Fig. [3] for axis orientations 
i. Open symbols are experimental measurements; error bars are comparable to or smaller than 
symbols. Solid lines are linear fits. Triangles show MSD slope of 1. 

and 2. If the particles in the trimer are identical, the threefold symmetry axis of the trimer 
ensures that Dt^i = Di^ and D^^i = -Dr,2 [37]. Thus the differences between these elements 
of the tensor imply that the particles in our trimer are not in fact identical. We performed 
HYDROSUB calculations to confirm that weakly breaking threefold symmetry results in dif- 
ferences between the in-plane elements of D. Our measured ratio Dr^i/Dr^2 = 1-03 ± 0.02 
corresponds to a 3% size difference between the spheres. This is consistent both with the 
particle manufacturer's certificate of analysis as well as with particle size differences deter- 
mined from holograms. This shows that the measurements are precise enough to detect the 
weak breaking of threefold rotational symmetry due to particle polydispersity. 

Our work demonstrates experimentally how both large and small differences in the sym- 
metry of small particles affect the diffusion tensor. The technique we use, holographic 
microscopy, can measure elements of the diffusion tensor to high precision, 1% or better, 
small enough to resolve weak symmetry breaking due to particle polydispersity. The high 
precision is enabled by the inherently short acquisition times of the technique, which allow 
us to observe rapidly diffusing systems on time scales ranging from a fraction of a rotational 
diffusion time 1 / D^^i to hundreds. Although here we have measured the diffusion of isolated 
clusters, it should be possible to use the same technique to measure diffusion tensors in 
environments that are relevant to self-assembly but challenging for computation. In partic- 
ular, it may be possible to measure diffusion tensors near boundaries or other particles, or 
for clusters that have internal degrees of freedom. There are only a few studies examining 
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TABLE II. Measured diffusion tensor elements for trimer shown in Fies. [3] and ID with comparisons 
to computations from HYDROSUB [32j. Computations use a = 500 nm obtained optically from the 
best-fit hologram models and rj = 1.049 mPas from single-particle diffusion data; the difference in 
r) from the dimer measurements is due to a difference in room temperature. 
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Dr,2 (s-1) 


0.270 ± 0.002 


0.296 


Dr,3 (s-1) 


0.210 ± 0.002 


0.220 


Dr,l/Dr,3 


1.32 ± 0.02 


1.34 


Am/A,2 


1.03 ± 0.02 


1.00 


A,i (x 10-13 m^s-i) 


2.466 ± 0.015 


2.64 


A,2 (x 10-13 m^s-i) 


2.446 ± 0.015 


2.64 


A,3 (x 10-13 m^s-i) 


2.372 ± 0.015 


2.41 


A,i/A,3 


1.04 ± 0.01 


1.09 



the effect of interparticle hydrodynamic couplings on particle diffusion, and these have been 
restricted to spheres in planar geometries [381 139] . Furthermore, measurements on other, 
even less symmetric clusters may be able to reveal translational-rotational coupling as well 
as the off-diagonal rotational and translational elements in the diffusion tensor. 
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W. Perry for helpful discussions. This work was supported by the National Science Founda- 
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Supplementary Information 

Hologram Analysis and Data Reduction 

We describe in further detail how we fit scattering models to the holograms we record of 
dimers and trimers to extract 3D dynamical information, with attention to issues that arise 
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with the large number of holograms (~20,000) we must analyze for each case. 
Model Fitting Procedure 

Our technique of fitting scattering solutions to holograms has previously been described 
fl9\ |20] . For the dimers, we fit a scattering model that depends on one refractive index for 
both particles, the radius of each particle, the 3D position of the dimer center of mass, 2 
orientational Euler angles, and a scaling parameter (TH]. The model for trimers differs only 
in that we fit for three orientation angles and for only one average radius [20] . 

The largest bottleneck in fitting scattering models to large numbers of holograms using 
Levenberg-Marquardt minimizers is that an initial guess for the model parameters needs 
to be provided for each hologram. With time series of holograms, only one initial guess is 
necessary, in principle: we can use the best-fit parameters of each hologram as the initial 
guess for the next. We have found it effective to speed up this process by doing an initial 
rough fit to a randomly chosen subset of 10% of the pixels of each hologram. Subsequently, 
we use the parameters obtained from the rough fits as initial guesses for a fit to all the pixels. 
We do this second stage of fitting in parallel. 

Validating Model Fits to Trimer Holograms 

The model for the trimer holograms has an additional orientational degree of freedom 
compared to the model for the dimer holograms. We have noticed that on occasion the fit- 
ter converges to best-fit parameters that result in the best-fit model hologram having subtle 
differences when compared to the experimental hologram; this usually stems from the orien- 
tation angles being incorrect. We do not observe this problem for the dimer holograms. To 
detect holograms with potentially incorrect best-fit parameters, we inspect the statistic 
[20] of the fits. We also compute a statistic for a binary version of the experimental 
and best-fit holograms, where all pixels above the mean of 1 are set to a value of 1 and all 
remaining pixels are set to a value of 0. The binary image is much more sensitive to the 
shape of the interference fringes. 

When we compute correlation functions such as mean-squared displacements from the 
trimer holograms, we reject the contribution from any holograms where either or binary 
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is worse than 2 standard deviations from a rolhng mean. Manual inspection of 200 
randomly chosen trimer holograms that were not rejected under these criteria revealed 7 
questionable fits. We infer from the Poisson distribution that, to a 99% confidence level, 
the percentage of remaining bad fits is less than 8%. We also reject the contribution from a 
given pair of holograms if the probability of obtaining either a center of mass displacement 
or angular displacement of the observed magnitude is less than 10~^. We compute these 
probabilities using estimates for the diffusion tensor elements, and choose the threshold of 
10^^ to avoid biasing the observed distribution and to make the cutoffs weakly sensitive to 
the estimates for D. 

Performing this cutoff procedure requires knowing the probability distributions governing 
translational and rotational displacements. The probability distribution for translational 
displacements is Gaussian, but the distribution function for rotational displacements is not. 
Instead, the probability density function fi{9; r) for observing an angle 9 between cluster 
axis Uj(t) at a given time t and after a time interval r is given by 



where -Dr-,eff = {Drj + Dj.^k)/'^-, Drj and Dr^k are the elements of D^'' describing rotations 
about the two cluster axes other than i, and Y^{6) denote spherical harmonics with m = 0. 



As a final verification that our holographic imaging is correct and that any remaining 
errors do not substantially affect the dynamics we measure, we compute probability distri- 
bution functions for the dynamical quantities we use to measure D from the data. Figure 



find that the measured distributions agree well with the expected distribution computed 




oo 




We briefly discuss the origin of this distribution function in Section [} 





Estimating Solvent Viscosities 



Here we discuss in greater detail the inference of the solvent viscosities rj needed to com- 
pare experimentally measured elements of the diffusion tensor D to theoretical predictions. 

13 




4 L U I ILl I I I I I ijj 

-3 -2-10 1 2 3 



FIG. 5. Distribution functions for trimer angular displacements and cluster- frame displacements. 
Histogram points computed from experimental data are shown in open symbols; solid lines show 
theoretical predictions computed from elements of D reported in Table II of the body of the paper, 
(a) Rotational dynamics of U3. Predicted distribution computed from Eq. [5] (b) Cluster- frame 
displacements along axis 3. Theoretical distribution is a Gaussian with a mean of and a variance 
of 2A,3r. 



The viscosity of our solvent, a mixture of H2O and D2O chosen to density-match 
polystyrene particles, has a strong temperature dependence. Figure [6] shows that the 
viscosity of a bulk sample of the solvent, measured using a Cannon-Manning capillary 
viscometer, varies by nearly 20% over a 6°C temperature range. We performed all the 
experiments described in the body of the paper at room temperature, but we have observed 
that the room temperature can change by several °C over the course of a few hours, most 
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FIG. 6. Temperature dependence of solvent viscosity. Data points, open symbols, were measured 
with a Cannon-Manning capillary viscometer. The solid line is a best-fit quadratic function that 
allows for interpolation between the measured points. 

likely due to the cycling of the building heating and air conditioning systems. Moreover, 
particularly if the laboratory room temperature is changing, the temperature in the sample, 
sealed in a glass sample cell, may differ from that of the surrounding air. 

Consequently, we believe that that the best way to estimate the solvent viscosity is to 
observe the in situ diffusion of single colloidal spheres, which are always present in the sample 
due to the arrested aggregation technique we use to make the clusters, either immediately 
before or immediately after imaging the diffusion of a cluster of interest. The Stokes-Einstein 
relation gives the translational diffusion constant |10] D in terms of the temperature T, the 
particle radius a, and the solvent viscosity rj: 

6nr]a ^ ^ 

Using Eq. [6} once we determine D for a diffusing sphere of radius a, we can infer the ratio 
ksT/r]. Because of the strong temperature dependence illustrated in Figure |6| ksT and r] 
should not be viewed as independent parameters. Moreover, from dimensional considera- 
tions, the elements of D are always proportional to ksT /rj. Once we determine ksT /r], we 
use the best-fit line to the data in Figure [6] to infer rj and ksT separately. While this is not 
the usual context in which microrheological experiments are performed, we essentially treat 
the diffusing single spheres as in situ thermometers. 

We obtain D from an MSD computed from the 3D trajectory of a diffusing particle: 
(Ar^(r)) = QDt. In all cases, we obtain the trajectory using DHM and record holograms at 
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25 frames per second. We obtain a radius, index of refraction, and 3D position from each 
hologram by fitting a model based on the Lorenz-Mie solution [TH] . 

For the dimer experiment, which used particles with a nominal radius of 650 nm, we 
measure D = 2.533 ± 0.017 x 10^^'^ m^s^^ for a diffusing particle with an optical radius of 
639 nm. If we assume that the particle has the same enhanced hydrodynamic radius of 709 
nm as we inferred from the dimer data, independent of any considerations of ksT or rj, we 
can subsequently use the data in Figure [6] to infer a solvent viscosity of 1.187 mPas, which 
is within 3% of the best-fit solvent viscosity, 1.159 mPas. The consistency of these values, 
along with the excellent agreement between the measured and predicted values of D\\/D±, 
which is independent of a and kBT/rj, validates our dimer measurements. 

For the trimer experiment, we measured D = 3.996 ± 0.055 x 10~^^ m^s~^ for a diffusing 
sphere of nominal radius 500 nm. With no analytical theory as we had for dimers, we cannot 
rigorously find a best-fit radius for the trimer. We take the optical radius of the particle, 
517 nm, as an estimate of the particle size and use the data in Figure [6] to infer 77 = 1.049 
mPas, the value we use in the hydrosub calculations. 

Probability Distribution for Rotational Displacements 

We briefiy describe the origin of Eq. [5} the probability density function for finite rotational 
displacements. As described in Section [} we use this distribution to reject pairs of holograms 
that exhibit highly improbable angular displacements, most likely due to an incorrect model 



To calculate the relevant distribution function, we will consider what happens to an imag- 
inary ensemble of clusters undergoing rotational diffusion. We will assume that translation- 
rotation coupling is negligible, so that we can consider the rotational motions independently 
of the translations. Suppose that we observe only the motion of one body the 
clusters in the ensemble undergo rotational diffusion. Lastly, suppose that that we prepare 
the ensemble such that at t = 0, Uj for every cluster lies at the same point on the unit 
sphere, which we may choose to be at = without loss of generality. We seek to compute 
the probability distribution fi{0,(f); r) such that 



fit. 




(7) 
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gives the probability of finding Uj between 0o and 0i and between Oq and at t = r. 

For isotropic rotational diffusion, such as that of a sphere, calculating /j is straightforward: 
the diffusion can be described by a rotational Pick's law characterized by a single rotational 
diffusion constant Dr [S]. The initial condition 



where S{9) denotes the Dirac delta function, then determines fi{0,(f); r). This idea can be 
generalized to the case we are interested in, where D'"'" is diagonal but not isotropic, using 
Brenner's tensorial formalism [22j. Because the details are quite involved j42j . we instead 
give a physical argument that allows us to apply the isotropic solution to the anisotropic 
case. 

The time evolution of must be governed by Drj and -D^. ^, the diffusion constants for 
rotations about the other two cluster axes. In our ensemble of clusters, prepared such that 
all clusters initially have Uj at ^ = 0, the clusters will not all have the same orientation: 
and Ufc can lie anywhere on the equator of the unit sphere. Consequently, observing only 
the motion of Uj, we will on average observe evolving according to an effective rotational 
diffusion constant -Dr.cfr where -Dr,efT = {Drj + Dj.^k)/'^- We may then straightforwardly adopt 
the result from isotropic diffusion [41j, which leads to Eq. [s] Note that fi is ultimately 
independent of because of the symmetric manner in which we prepared the ensemble. 




(8) 
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